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I. INTRODUCTION 


In order to use perturbative QCD (pQCD) for the description of exclusive processes, one 
usually appeals to factorization theorems, which ensure that the pQCD calculable infor¬ 
mation, described by the so-called hard scattering amplitude, Th, can be factored out at a 
factorization scale /i|, whereas all pQCD non-calculable (i.e., non-perturbative) contributions 
are encoded in terms of universal hadron distribution amplitudes (DA)s, tfh- Schematically, 
this can be illustrated for the case of the 7 * 7 * —>■ vr transition form factor in the following 
way: 


q^) = UTh{x; /ip) ® l4) > 

X 


( 1 . 1 ) 


where and are the photon virtualities and 0 denotes the usual convolution symbol 
{A{x)®B{x) = L dxA(x)B(x)) over the longitudinal momentum fraction variable x. The 

X 

dependence of the DA on the factorization scale is governed by the Efremov- 

Radyushkin-Brodsky-Lepage (ERBL) evolution equation [l|, 


d In /i| 


V(x, u; 0 /ip) 

U 


( 1 . 2 ) 


with the evolution kernel 

V{x,u■as{^l^))) = V2{x,y) + ..., (1.3) 

adopting the notation of j^. The one-loop evolution kernel Vq was introduced in Ref. j^; 
an analogous expression for V 2 at the two-loop level was derived in Refs. 0, S 0 01- 
factorization applies, it is, in principle, safe and legitimate to use different values of /i|, 
dispensing with the factorization-scale dependence by means of the renormalization group. 
Usually, one sets yp = in order to eliminate in Th large logarithms of the 

type ln((5^//i|) and for that reason one then needs to evolve the meson DA in 

accordance with Eq. (O from a low normalization point /Xq ^ 1 GeV^, at which the pion 
DA has been (non-perturbatively) determined, to a higher scale y? = at which comparison 
with experimental data may be attempted. 

The structure of the two-loop evolution kernel was analyzed by Mikhailov and Radyushkin 
(MR) in ^ and its diagonality violating terms (in the Gegenbauer basis) were identihed, 
though their origin remained partly unclear. Moreover, the two-loop evolution with of 
any pion DA at the reference (or initial) momentum /Xq was carried out numerically, but 
including only the hrst few expansion coefficients at order a^. On that basis, these authors 
concluded that the next-to-leading-order (NLO) correction to the pion DA with /Xq ~ 1 GeV^ 
and at a; = 0.5 remains less than a few percent even for momenta far beyond 10^ GeV^. 
Later, the two-loop evolution of the pion DA was studied by Muller in using conformal 
constraints. In this work, a complete formal solution in NLO was obtained with the inclusion 
of all two-loop mixing coefficients. This solution was further discussed in for a running 
and also a hxed coupling constant. It was found that the NLO correction can be rather 
large - especially for endpoint-concentrated DAs, like the Ghernyak-Zhitnitsky (GZ) na 
one, supplying logarithmic enhancement exactly in this region. 

However, as we will show below, all these solutions violate the group character of the 
renormalization-group (RG) evolution transformations to the order Whether it is possible 
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to obtain an approximate solution of the evolution equation that, nevertheless, respects its 
group properties, remains an open question. It is exactly this issue to which the present work 
is devoted. Such an improvement is not only theoretically important, it is also of practical 
concern since the modihcation of the evolution behavior, entailed by the restoration of the RG 
properties, will have influence on measurable hadronic observables, like meson form factors. 

In the next section, we shall briefly review the main properties of the two-loop approxima¬ 
tion of the QCD coupling and describe the standard formalism for the DA evolution at the 
two-loop level. Section IIIII presents the analysis of the RG transformation when the second 
and higher Gegenbauer harmonics are taken into account. A semi-explicit solution of the 
two-loop evolution equation that preserves the RG structure is constructed in Sec. IIVI In 
Sec. El we discuss the numerical importance of the obtained solution. The special case of 
a fixed coupling constant is treated in Sec. ED while some important technical details are 
collected in two appendixes. 


II. QCD COUPLING AND EVOLUTION OF THE PION DA AT NLO 
A. Standard formalism with standard notations 


To illustrate the “RG philosophy” and clarify the objective of this work, let us recall the 
Ovsyannikov-Gallan-Symanzik equation for the (running) coupling, as (Q^), in QGD: 


das (Q^) 
dln(g2) 


(d{as{Q‘^)) . 


( 2 . 1 ) 


Strictly speaking, asiQ"^) and also /3(as(g^)) depend on the number of active flavors, Nf. For 
the considerations to follow, this will not be important and hence we omit the ^/-dependence 
any further. The /3-function in the NLO approximation is given by 

' 3 (“*) = -s(‘» + '’‘S)' P' 2 ) 

where the standard /3-function coefficients are provided in Appendix El The two-loop equa¬ 
tion for asiQ"^) (with Aqcd = A) reads 


dvr 


boas{Q‘^ 


Cl In 


dTT 


hoasiQ"^ 


+ Cl 


= In 


A 2 


The solution of this equation can be written in terms of the 

l/ci 


argument ((Q^) = 


to obtain 


with Cl = 61/60 • (2.3) 

Lambert IF_i-function with the 


a 


2-loop 


(«6 = 


dvr 


i>o(iV/)ci [l + M-'-iKW"))]’ 


(2.4) 


as it was shown in (TT. . 

Now let us turn to the ERBL evolution equation for the meson DA. First, we recall the 
main properties of the one-loop approximation. (See Ref. Q for a pedagogical exposition 
and further references.) The one-loop evolution kernel has a very simple structure, viz.. 




Vi.ioop{x,y]as) = ■^Vo{x,y) 


(2.5) 
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with factorizing Og- and x-dependences. The formal solution of the ERBL equation in the 
one-loop approximation is 


= exp 






47r/3i(a 


-da ® 


( 2 - 6 ) 


where Pi{a) = —boa‘^/{4:7i) is the one-loop /3-function and the exponent above has to be 

r 1 ^ 

evaluated according to Vq ® ipn{x; = Vo{x,ui) <S) ■ ■ ■ Vo{un-i,Un) ® It 

L J lip Ufi — i Uji 

is useful to expand the meson DA^ 

‘P7,{x;Q^) = Q{x)'^ aniQ'^) - tpn{x) 


(2.7) 


n>0 


in terms of the eigenfunctions Q{x)ijjn{x) of the one-loop ERBL kernel i.e., in terms of 

the Gegenbauer polynomials 

f2(x) = 6x(l —x), ipn{x) = G^^^(2x —1), (2.8) 


corresponding to the eigenvalues 

^ (^) 7o(»j) (2.9) 

with 7 o(n) given in Appendix El In this representation all the dependence on is contained 
in the coefficients a„(Q^): 


l-loop 


(Q^) 



' as{Q^) 


7(n) 


with 7 (n) 


7o(n) 

2bo 


( 2 . 10 ) 


This simple scheme is due to the special factorized structure of the one-loop evolution kernel 
()2.5|) . In the two-loop approximation this is no more the case. Indeed, the evolution kernel 

I4-ioop((r,?/;a,(/i^)) = Voix,y) + V2{x,y) (2.11) 


has, at each scale /i^, different eigenfunctions, which explicitly depend on Q;s(/i^). Never¬ 
theless, these eigenfunctions of the two-loop ERBL kernel can also be expanded in terms 
of the one-loop eigenfunctions r2(x)'0„(x). In this basis, the two-loop kernel ( 12.1111 can be 
represented in a matrix form of a triangular type: 


h2-loop((^) y 1 ^s) 


''2-loop 


(tts) 


^»((g) ( 2 . 12 ) 

n j ^ 

~Y (S^) + (^) “ ^hn0{j < ri)]} (2.13) 


^ Here denotes the sum over even indices n > 0 only in order to account for the symmetry relation 

n>0 

(pTr(x; Q^) = ipTr{l — x) due to charge-conjugation invariance and isospin symmetry. 
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with normalization coefficients Nj, next-to-leading order anomalous dimensions 7 i(n), and 
off-diagonal matrix elements Mj^n given in Appendix ^ . 

We see that the diagonal terms in Eq. (EH, ^2- 'ioop('^s), define the anomalous dimensions 




-1 



loin) + 



(2.14) 


alMj^n, define the mixing of higher harmonics. 


whereas the off-diagonal terms, 

Solutions of the evolution equation at the two-loop level have been given in Eiia , and 
have the form 


n 


Mx) + dn,j{Q\ ^4)i^j{x) 


j>n 


(2.16) 


SO that the evolved coefficients are 


(2.16) 


0<j<n 


Here, the “diagonal” part E„(Q^,/iQ) is the exact part of this solution, namely. 


EniQ^d) = 

Si = 


en(Q^ 

hi 

Auhr) ’ 


26 ’ 


r 1 7 (n) r 

^niQ'^) = oisiQ"^) l + Siasi^Q"^) 

, X 7iW^o - 7o(^)&i 
-’ 

2bobi 


uj{n) 


(2.17) 

(2.18) 


while the “non-diagonal” part is considered in the NLO approximation. Note that the coef¬ 
ficients dnjiQ'^, fXo) are related to the off-diagonal matrix elements (I2.1dj) and fix the mixing 
of the higher, j > n, harmonics in 

In the MR solution these coefficients are given by 




M,,„ j ^ 

26o[7(?i)-7(i) - 1] 1 L«8 (/^o). 


7W-7(i)-l' 


(2.19) 


As we shall see in the next sections, the approximate MR coefficients (Q^, /Iq) explicitly 
violate the group property of the RG transformation at the 0(Q(s)-level. 


B. Standard formalism with modified notations 

It turns out to be more convenient to rewrite the main equations of the previous subsection 
in terms of the modified notations, ■j{n) (Eq. (I2.1()j] h di, and uj{n) (Eq. ()2.18j) b Then, the 
/3-function becomes 

fS (os) = ~ho—^ [1 -|- diQ;<j] . (2.20) 

4:71 

Next, the two-loop kernel reads 

KiLpM = -bo I [7(n) (1 + hia,) + u{n)6ias] Sn,j - < n)| . (2.21) 
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Finally, we rewrite the two-loop ERBL equations (Cl, (EH, using the representation 

(d and changing variables from to as with the help of Eqs. m to derive 


= -Y (S) 'Mj,na,{as) 


0<j<n 


( 2 . 22 ) 


After taking into account Eqs. TTM and dm, this expression gives 

as (1 = [-)/(n) (1 u!{n)6ias] an(a^) - ^ 7hn,jaj{as ), (2.23) 


0<j<n 


where we have defined the reduced matrix elements runj = Mj^n/{87^b(j). These expressions 
are going to be helpful in deriving an exact RG solution in Sec. IIVI 


III. GROUP PROPERTY OF THE RENORMALIZATION GROUP TRANSFOR¬ 
MATION FOR GEGENBAUER HARMONICS 

In this section we consider exclusively the two-loop evolution equation given by expressions 
(d and (cnD- For this reason, we omit in the following all two-loop-superscripts. We 
want to understand if the approximate form respects the group property of the RG 

transformation, that is. 




(3.1) 


if not, we want to estimate to what extent it violates it. 


A. The second harmonic 

To this end, let us consider first the evolution of the Gegenbauer coefficient a 2 - According 
to ()2.1f)|l . we have 

a2(/^o) ^ ®2(Q^) = t^(Q^,ho)®2(/io) = -S2(Q^,/^o) ®2(/^o) +-^2 o(Q^, ho), (3-2) 

where we defined D 2 o{Q‘^,Y) = (Q<s(Q^)/4vr)(io2(Q^, h^)- If group property (j3.1|l is valid, 
then it follows 

E2(Q^ q^) [E2Y, h')a2(/i') + D20Y, h')] + ^2o(Q^ q") = ^2(Q^ h')a2(h') + /^2o(Q^ h') 


and from the arbitrariness of 02 (h^) one gets 

E2{Q^Y) = E2{Q\q^)E2{q\Y)-. (3-3) 

D2o{Q\Y) = E2{Q\q^)D2o{q\Y) + D2o{Q\q^). (3.4) 

The first equation is satisfied identically by virtue of Eq. (cnD. The second one is more 
complicated, but can be readily proved. Define a function ^ 2 o(Q^h^) by 

I^2o(Q^,h^) = G2{Q‘^)Z2o{Q‘^ , Y) ■ (3.5) 

Then, Eq. (13.4|) implies 

Z2o{Q\Y) = Z2o{Q\q^) + Z2o{q\Y)- (3-6) 
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The general solution of this functional equation is just 

Z2o(g^/x^) = T2o(g^) - T2o(/i^) (3.7) 

for some auxiliary function T 2 o(/^^). In order to hnd the explicit form of \h 2 o(/i^), one has to 
use evolution equation (EH, what we will conduct in Sec. nYi By using the MR approximate 
solution (EIHD, we see that it generates the function 

OjTWtf*") = -'•' 20 (^ 1 ")] (3-8) 

with 


\jf f ‘ 2 \ _ ^^^20 

2«('* ) - (^(2) - 1) |a.(^2)J,(2)-l ■ 


(3.9) 


From this, we conclude that in order that is consistent with the RG transfor¬ 

mation property, it has to be given by 


= e^m 'l'2o(0") -'l'2o(/i") 


(3.10) 


Comparing (EHl) with (EUB leads us to the conclusion that the group property of the RG 
transformations in the MR approximate form is violated because the factor [l-|-5iQ!s(g^)]^*'^^ 
is missing. Evidently, this violation is of 0{as). 


B. The reth harmonic 


The next task is to generalize these results to the case of arbitrary polynomial order 
n = 2,4,.... To achieve this, we define Dnk{Q‘^, and ZnkiQ"^, as 

= ^4^£,(<3^^^")d,:„(<3^^^") = e„(<3=)Z„.j(Q2,P)e,(;i2)-u (3,11) 

= ^^ei(<2")4»Wtl‘te„(<2")-‘. (3.12) 

47r 

Then, the evolution of the Gegenbauer coefficient from the scale to the scale in 
accordance with Eq. is given by 


an(g^) = e„(g^) 


(3.13) 


0<j<n 

The group property EH) dictates 

[^r^J(g^g^)+ ^r^J(g^/i^)] ai(/i^)e“^(/r^) 

0<j<n 

+ E ’ Z Tj(<2t?")2?,i=(?tl<tat(l‘tei‘(l‘") (3.14) 


0<j<n 


0<j<n 0<k<j 

and from the arbitrariness of we obtain 

■^n,j(g^,h^) = Znj{Q‘^,q‘^) + Znjiq'^, ^‘^) + ^ ZnkiQ"^, Zkj{q^, ■ (3.15) 

n>k>j 
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1. The case j = n — 2 


For j = n — 2 we have the complete analogue of Eq. (ESD, i.e., 
Zn,n-2{Q‘^ 1 = Zn,n-2{Q‘^ , (f) + Zn,n-2{(f , 

yielding the exact solution 

Zn,n-2iQ^^ = ll'n,n-2(Q^) ~ d^n,n-2(h^) • 


(3.16) 

(3.17) 


2. The case j = n — k with arbitrary and even k < n 


Let us rewrite Eq. (ITT3D in the more appropriate form 

Zn,n-k{Q‘^, = ^n,r^-fc(Q^g^) + ^n,n-fc(g^h^) 

+ ^ Zn,n-j{Q‘^,q^)Zn-j,n-k{q^,h?) ■ (3.18) 

0<j<k 

Using the principle of mathematical induction, we can prove that 

Zn,n-k{Q^,T^) = '^n,n-k{Q^) “ d^n,n-fc(/i^) “ ' Zn,n-j{Q^, n-j,n-k{T^) (3.19) 

0<j<k 

is the solution of Eq. (I3.18|l (relegating further details to Appendix lB|l . 

We are now in the position to summarize our Endings and rewrite the solution of the NLO 
evolution equation in the form (EH) with the nth Gegenbauer coefficient that respects 
the RG properties being given by 




an(/io)en^(ho) + (F^/^o)«i(ho)e/(ho) 

0<j<n 


where 


r "17(’^) r "I 

e„(/i^) = a^(/i^) l + 6ias{T^)^ 

ZnSiT^Tl) = ^n,fc(h^) - ^n,fc(ho) - Y1 /^o) (ho) . 


1 


n>j>k 


(3.20) 


(3.21) 

(3.22) 


with 'h„j(/r^) being defined on account of evolution equation ()2.23jl . a task we will conduct 
in Sec. IIVI 

It is instructive to rewrite the approximate MR solution in an analogous manner to obtain 


ar(h^) 


en(h^) 


Q'n(P'o)®n (ho) A ^ ^ ^n,j (h )ho)®i(ho)G’ (ho) 


0<j<n 


efc(/i^)Qi,(/i^)^(^^ r^MR/ 2 n _ ^MR/ 2^1 

e„(/i2)«,(/i2)7(A:) \-^n,k{T ) ^n,fc (ho)J 

[1 + rah^)-^»h^)] , 

'Tklrijn—k 

|^(„) _ ^(„ - i) - 1] |a.((32)]’<”l-''(”-'=>-‘ ■ 


(3.23) 


(3.24) 

(3.25) 


Gomparing Eqs. ()3.22j) and (13.2411 . we see that Z^^(/i^,pQ) differs from /Iq) in two 

respects: 
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• It contains a factor [1 + in the leading term — d'n,fc(/^o)- 

• It comprises additional terms of the type pQ)\I'j^A:(/^o )5 which are related to j 

bearing a smaller index j < k. 


IV. EXACT RENORMALIZATION GROUP SOLUTION 


Up to this point we have presented only an approximate form of the d'-functions (see, 
Eq. ()d.25|l h derived from the MR solution, that satisfy the RG equation. Now that we have 
outlined the basic steps of the RG restoration, we proceed to our analysis of its consequences 
for the \I'„_A;(/i^)-functions. Substituting our solution (Id.2211 into the two-loop ERBL equation 
(IT^ . we obtain^ 


(1 -|- ^ ^ 


dZnj (as,ao) 


' Mj,n 


0<j<n 


Sirbc 


0<j<n 

ej{as) 


da. 


aj{ao)e- (oo) — 


aj{ao)ej^{ao) + ^ Z'^^{as, ao) ak{ao)e^^{ao) 


0<k<j 


(4.1) 


Ghanging indices of summation in the double sum from j, k to k,j and rearranging them as 
to sum Erst over j, we obtain due to the arbitrariness of aj{ao): 

dZ^ j (o^s? ^o) ... ^^ ^ RO 

(1 -h 5 iQ!s) en{as) --= -mn,j ej{as) - ^ mn,k («s, ao)ek{as) . (4.2) 


For j = n — 2 this differential equation reduces to 


j <k<n 


dZ'^^_2{a,ao) 

da 


whose exact solution is 


^Sn-2(«,ao) = ^n,n-2(«) “ ^n,n-2(ao) 


^ri.,n —2 


fRG 


en- 2 («) 


en{a) (1 6ia) ’ 


RG 


with 


2 (a) ; 


,RG 

,n 

sRG 


(4.3) 


(4.4) 


(4.5a) 


crhg ^ - 2^1 (1 - 7 (n) + 7 (j), 1 + uj{n) - u{j), 2 - 7 (n) + 7 (j), -6ia) 

By definition, the function <l>J^^(a) satisfies the following evolution equation: 


da 


m. 


i,j g('^) 


enia) (1 -h ^iQ?) 


(4.6) 


In Fig. [T] we show both solutions, Zf-Q^a, Gq) and Z^^{a, ao); for a < ao = 0.5. We observe 
that the error in Z 2 fl{a,aQ), when using the approximate MR solution, varies from 25% 
(a ~ ao) to 15% (a —>• 0.1). 


^ From now on we trade all the /i^-dependence for an ag-dependence (with ao = as(/j,g)). 
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0.1 0.2 0.3 0.4 0.5 

a 


1 

0.95 
0 .9 
0.85 
0 . 8 
0.75 
0.7 

0.1 0.2 0.3 0.4 0.5 

a 



FIG. 1: (a) The solid line corresponds to (a, oq) (RG-improved case) and the dashed line to oq) 

(MR-approximate case), (b) The ratio Z^^{a,aQ)/Z^^{a,ao) is plotted vs. a. In both panels we use 
oo = 0.5. 


Consider now Eq. ( 1121 ) in the general case k = n — j with 2 < j < n. We have the RG 
representation of its solntion dm, which can be snbstitnted into the evolution equation 
dOl) to obtain^ 


da 


= —m 


ej{a) 


E' 


”’^e„(a) (1 + (5ia) (1 + ’ 


^fc,j(a)efc(a) 


(4.7) 


The solution of this equation can be represented in terms of quadratures: 

i"'n,k I - T 

L e„ I ( 

j<k<n 


Jo 


d a. 


/o e„(a)(l + hia) 
Before analyzing this solution, let us sketch the procedure to derive it. 


(4.8) 


• We analyzed the NLO evolution equations 

• We used for the representation (Id. 2011 in terms of pg) and obtained Eq. 

(I121)- 

• We employed the RG representation for and derived this way Eq. 

(EZZD which yields for an explicit solution given by expression (USD- 

An exact solution of this equation for the case fc = n — 2 is provided by 
(I4.5bfl . expressed in terms of the hypergeometric function 2 F 1 {a,b,l + a, —6ia) with a = 
1 — 7 (n) + 7 (n — 2) and 6 = 1+ a;(n) — a;(n — 2). From Eq. (I4.8jl we can restore the leading 
asymptotics of 4 /^j(q;) to get 



a^O + 


(4.9) 


a form dictated by the <I>^^(Q:)-term, whereas the term containing the Yl' Eq. (14.8 |) gen¬ 
erates a correction with the asymptotics _ TPjg property allows us to suggest an 


^ We use here the same trick of interchanging summation indices (j, k) —> (fc, j) and performing the sums in 
j and k, as in Appendix IBl 
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1.3 

dS/dS 

1.3 



. 



1.2 

■ ■ ■ 

1.2 


1 .1 

■ 

1.1 

rjRG.app / 7-)RG 
’ ^n,0 / ^nfl 


;-)RG,app / piKG 
^n,0 / ^n,0 



t * 4 ♦ ♦ ♦ ♦ ♦ 


* ♦♦♦♦♦♦♦ 






± ^^^^^^_ 




2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 20 

n n 


FIG. 2: Comparison of the RG improved, RG approximated, and original MR evolution in terms of the 
coefficients I?„^o(ai, cto) with ai = as(4 GeV^) and ao = as(l GeV^). The blue points at the bottom denote 
ao), whereas the red points above them represent '^o)- Left 

panel: The original values of Mk^n with Mo ,4 = 0.285 are used. Right panel: All coefficients Mk^n have the 
original values, except for Mo ,4 for which we set Mo ,4 = —0.8. 


approximate RG solution of the NLO evolution equation (indicated below by the superscript 
‘app’); viz., 

« 0 ) = . (4.11) 

The results for the coefficients Dnfi{ai,ao) are shown in Fig. |21 As already mentioned, the 
RG-approximate evolution is rather good for n > 6, but for n = 4 it drops to a much too low 
value compared to the exact RG-improved result. The reason for this “jump” can be traced 
to the (unexpected) smallness and positiveness of Mo ,4 = 0.285, relative to Mq ^2 = —6.01 and 
M 2,4 = —17.05, cf. Eq. (lA.Of) . To show that this jump is connected to the particular value 
of Mo,n = 0.285 at n = 4, and is not a numerical error, we show in the right panel of Fig. |2l 
a graphics, where we set by hand the value of Mo,n = —0.8. As one sees, the result of this 
rough simulation yields to a complete cancellation of the dip at n = 4 turning it into a slight 
bump at n = 6, proving the consistency of our numerical algorithm. 

The numerical solution of Eq. (USD proceeds through the following step-by-step procedure: 

• Given that we know the solution for j = n — 2—supplied by 4 /J)'*)^_ 2 (q;) = <F^^_ 2 (a), cf. 

Eq. (I4.5bll — we determine n- 4 (®) solve numerically Eq. (gSl) for j = n — 4. 

This yields 






'^n,n—2 


^n-2,n-4(a)en-2(a) 
en(a) (1-b hia) 


• We then solve numerically Eq. (gSl) for j = n — k using the results 

• • •: ^n- 2 ,n-fc(«)} oi the previous steps to get 

• Finally, we solve numerically Eq. (USD for j = 0 , employing {\h^^(Q;),..., 
to find 4/^®(a). 

This procedure is repeated for all values of n, starting with n = 2, continuing with n = 4, 
and so on. The upshot of this procedure, ao), is shown in Fig. 121 in comparison with 
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FIG. 3: (a) The solid line corresponds to Zf^^{a,ao) and the dashed line to Z^^{a,ao). (b) The ratio 
^ 4^0 (a,a!o)/^ 4 ^(]^(a!,a!o) vs. a is shown. We use in both panels ao = 0.5. 


the approximate MR solution cto), using in both cases a < cto = 0.5. One notes that 

the error in Z 4 ^o(a, ao) induced by the MR approximation is quite substantial, varying from 
+20% (a ^ ao) to -5% (a ^ 0). 


V. NUMERICAL IMPORTANCE OF THE RENORMALIZATION GROUP SOLU¬ 
TION 


Let us now outline the importance of the RG improvement in terms of the evolution of 
the one-loop asymptotic DA using the MR scheme and the RG-improved one in comparison 
and continue then to the more complicated case of a double-humped, endpoint-suppressed 
pion DA obtained in 1^. The asymptotic solution to the one-loop evolution equation reads 


= 1 GeV^) = (fias(x) = 6x(l - x ); an(ho) = 


if n = 0 
if n % 0 


(5.1) 


1.25 

1.2 

1.15 

1.1 

1.05 

1 

2.5 5 7.5 10 12.5 15 17.5 20 

n 

FIG. 4: Comparison of the RG-improved evolution coefficients oo) with ai = as(4 GeV^) and 

ao = as(l GeV^) contrasted to the MR-approximate version, (a) The blue points denote D^^{ai,ao)^ 
whereas the red points mark ZI^^(ai,ao) (note the logarithmic scale used for the ordinate.) (b) The ratio 
Dn${ai,ao)/D^)^{ai,ao) as a function of the order number n is shown. 
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Then, in accordance with (EH, we have 




1 if n = 0 

if 71 ^ 0 ' 


(5.2) 


Figure l^a) shows the comparison of the two evolution schemes for the coefficients 
-Dn,o(Q^i fio) = 4 GeV^ with n = 2,..., 20, whereas the ratio of these coefficients 

is displayed in Fig. Hb). 

Two observations are worth noting: 

(i) The absolute difference GeV^, 1 GeV^) — GeV^, 1 GeV^) for all n is small 

(being of order —Eq. flT2D i. 

(ii) The achieved improvement is high, reaching a level of reduction of the evolution coef¬ 
ficients of up to 20%. 


The next step of the analysis involves the comparison of the results of the RG and MR 
evolution approaches as applied to the pion DA itself—Fig. 0— taking into account the first 
100 nontrivial terms, meaning summing up n = 2,4,.. .200 higher Gegenbauer harmonics. 
In panel (a) we display the function ^p{x)/x in the vicinity of the endpoint a: = 0, whereas 
in panel (b) we show the DA around the middle point x = 0.5. The key observation here is 
that in both panels, (a) and (b), the RG-improved approach produces slightly smaller results. 
More precisely, the DA gradient for the DA value turns out to be smaller by 4% at the origin 
and by 0.1% at the middle point x = 0.5. 

In our recent papers with S. V. Mikhailov 0, dealing with the extraction of constraints 
on the Gegenbauer coefficients of the pion DA from the GLEO data on the 7*7 —>■ tt 
transition form factor we used the standard MR approach to evolve the pion 

D A I 15 II from the normalization scale /Xq — 1 GeV^ to the scale /igy = 5.76 GeV^, introduced 
in [l8| by Schmedding and Yakovlev (SY), as being relevant for the GLEO experiment. These 
results are compiled in Tabled One appreciates that the evolution improvement here due to 
the RG modihcation is of minor importance, reaching just the order of 1%. 

In Fig. Ha) we show the same sort of comparison for the truncated inverse moment: 


{x ^)n = 3^aj(Q^). 
j<n 


(5.3) 



FIG. 5: The left panel shows the extreme endpoint region close to a: = 0 to effect the changes entailed by 
evolution. The right panel shows in units of as/(27r) the relative NLO correction using the RG-corrected 
evolution, where is defined as — yas(a:)] /^s.s{x)- In both panels, the solid line 

corresponds to the RG-improved evolution, whereas the dashed one represents the MR evolution. 
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TABLE I: Results of the evolution of the Gegenbauer coefficients a? and 04 from the scale fiQ = 
1 GeV^ to the scale /Xgy = 5.76 GeV^ for the asymptotic, BMS j3l) and CZ pion DAs. 


DA models 

02 (Po) 

afHdv) 


04 (^ 0 ) 

af^idy) 

af^idy) 

As 

0 

-0.004 

-0.003 

0 

0 

0 

BMS 

0.204 

0.144 

0.145 

-0.144 

-0.093 

-0.092 

CZ 

0.56 

0.403 

0.403 

0 

-0.004 

-0.005 


We see that in the RG-improved evolution scheme, {x~^)n approaches the limiting value 
{x~^)oo, represented by the dashed line, more rapidly. 

Note in this context that in a previous analysis with K. Passek-Kumericki and W. Schroers 
of the pion form factor we have analyzed the NLO evolution of the pion DA in the MR 
evolution scheme, using the property of the two-loop evolution just described. More precisely, 
if one is interested only in the value of the inverse moment, as this is exactly the case for the 
factorized part of the pion form factor, then it is actually enough to use for the calculation the 
LO evolution (inni), ensu ring this way an accuracy at the 1% level. Indeed, to establish this 
property, we analyzed in [l9| numerically the convergence of the truncated moment {x~^)n 
up to n = 100. In the present RG-improved approach, this property can be established even 
at lower values of n. 

Let us step one level higher and consider now the BMS pion DA 

7’bms( 2:) = 6x(l-x) 1 + (nl) € 2 ^“^{2x - 1) + (nl) {2x - 1) (5.4) 


as the initial input of evolution. Then, in accordance with (Em, we have 

= E2{Q\ d) aTHd) + D2,oiQ^ d) 

= E,iQ^ d) aTHd) + d) + d) aTHd) ■ 

All higher coelRcients at the initial point — l GeV^ vanish and therefore we have 


(5.5) 

(5.6) 


''n>& 


(Q^) = D„,o(Q^ho) +^n,2(Q^ho)«2^®(ho) +^n.4(Q^ho)«4“®(ho)• (5-7) 


3.172 
3.17 
3.168 
3.166 
3.164 
3.162 
3.16 

0 50 100 150 200 0 50 100 150 200 

n n 

FIG. 6: The impact on the asymptotic pion DA and the BMS one (indicated by corresponding acronyms) 
of the RG-improved evolution (blue squares on the top) relative to the MR evolution (red diamonds) on the 
truncated inverse moment Eq. ifOl . as a function of the number n of the Gegenbauer harmonics 

included. The limiting value (dashed line) is also displayed. 
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FIG. 7: (a) The effect of the RG-improved evolution (blue dashed line) on the initial BMS pion DA (solid 
line) is illustrated, (b) We show in units of as/i^-w) the relative NLO correction using on the initial BMS 
pion DA the RG-improved evolution, the emphasis being placed on the endpoint regions a; = 0 and a; = 1. 
Notice that A(p 3 j;;[g(a;) is defined as [</5 bms(2^; Q^) ~ 7^bms(^)] IVas{x). For the sake of comparison with Fig. 
13 we normalize this expression to (/?as(a;). 


Then, in NLO, the BMS pion DA evolved from the scale /Iq = 1 GeV^ to the scale within 
the RG-improved scheme, is given by 


= 6x{l-x) 




n>l 


(5.8) 


whereas its counterpart in LO, 7 ’bms(^)> obtained with the diagonal part given by Eqs. 
()5.5|) -()5.6p. 

The comparison of the results of the RG-improved and the MR-approximate evolution 
approaches, as applied to the BMS pion DA, is shown in Fig. [Staking into account the hrst 
100 nontrivial terms (n = 2,4,... ,200). In panel (a) we display the evolution effect of the 
RG-improved approach for the BMS pion DA. For the reader’s convenience, we use in our 
analysis the same numerical values Nf = 3 and = 4 GeV^ as in P|, but employing a higher 
initial point of evolution /i^ = 1 GeV^. The illustration of the evolution effect, especially on 
the endpoints (x ^ 0, x —*• 1), is shown in Fig. EKb). We will return to this issue and study 
the influence of the RG-improvement in our discussion of the evolution with a hxed in the 
next section. 


VI. RG SOLUTION IN THE CASE OF A FIXED COUPLING CONSTANT 

We turn now to a comparison of our results with those obtained in P] for the evolution with 
a hxed restricting attention to the asymptotic pion DA, the BMS case being analogous. 
To this end, let us adopt our forrnulae to this case, using for the corresponding quantities 
a tilde: en(Q^), and First of all, we realize that we have a different 

“diagonal” evolution 


rjiji) 


en(Q^) . 
en(/i^) ’ 
1 

2 Utt 


en(Q^) = 

+ (S) 




LQ2J 
7i(^) , 


( 6 . 1 ) 

( 6 . 2 ) 
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where Qq is an arbitrary auxiliary scale. Equations (EH), (IT^ . and remain valid, 

but the ERBL equation for an{Q^) gets modihed to assume the following form 

(6.3) 

0<j<n 


dan(Q^) 

dln(Q2) 


= -T]{n)an{Q^) + ^ 


Equation dOI) then transforms to 
dZ^(Q\l^^) 1 /o.,Y 


(ilnQ2 2 VdTT 

and for j = n — 2 it reduces to 


M. 




+ 


ek{Q‘^) 




j<k<n 


en-2{Q^) 

-.(-I 


dlng2 


2 Vdvr 


^n—2,n 


Its exact solution reads 


with 




en-2(g^) 


1 / Cr.s \ ^ Ln,n—2 


2 VdTT/ r]{n) - r]{n - 2) en{Q‘^) 

Ln.n-2 


a.. 


dvr/ 7 o(n) - 7 o(n - 2 ) + [^^/(dTr)] (71 (n) - 71 (n - 2 )) ’ 


d-‘n,n—2 Mn—2,n • 


(6.4) 


(6.5) 


( 6 . 6 ) 


(6.7a) 


(6.7b) 

(6.7c) 


But for this solution we have already obtained a RG-improved expression, given by dsm. 
After substituting this expression into the evolution equation (B, we obtain 


d^n,,{Q^) 

dlnQ2 


1 

2 



M. 


J,n 


en(g^) 


j<k<n 


efc(g^) 

en(g^) 


The solution of this differential equation can be represented as 




1 




e,m 


2 KAttJ r]{n) - r]{j) e„(g 2 ) 

lQ2fQ2f’‘>-nU) 


CXc 




An/ 7o(n) - 7o(j) + [a./( 47 r)] (71 (n) - 71 (j)) ’ 
where the matrix Lnj is dehned iteratively by 




1 


' 2 Vdvr; 


E /<> 


' L 

k.n 


k,j 


j<k<n 


Vik) - VU) 


( 6 . 8 ) 


(6.9a) 

(6.9b) 


(6.9c) 


In Fig.IHlwe show numerical results for Lnj and Mj^n- As one sees, the RG-improved solution 
Lnj is for 1-Cj<n — 23to4 times smaller compared to Mj^n- 
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FIG. 8: Comparison of Lnj and for fixed as = 0.5. (a) The case of a fixed difference n — j = 20 is 
shown, (b) An analogous situation is shown for n = 200. 


Having obtained these exact solutions, it is instructive to see how the result derived in 
for the asymptotic oo) solution of the NLO ERBL equation with hxed can be repro¬ 
duced. From Eq. we know that the Q^-dependence of the Gegenbauer expansion coef- 

hcients a„ is determined by the combinations e„((5^)e“^(/i^) and 
The hrst term has an evident asymptotics: 


en{Q^) 


q2^oo r 1, if n = 0 
to, if n 7 ^ 0 


( 6 . 10 ) 


whereas the asymptotics of the second term can be determined as follows. First, let us 
determine the asymptotics of this term in the case j = j 2 = n — 2: 




1 

2 KAtt) 




2 \ n{j2) 


For Q'^ —>■ oo we have 






1 / cuk \ ^ 


2 \ Vin)' 


^n,j2 


7]{n) - //(is) 


9 V /I J r V -^2 = 0 
2 KAtiJ r}[n) 

0, if j2 7^ 0 


( 6 . 11 ) 


The case j = = n — A differs from the case j = j 2 by the contribution of an additional 

term; viz.. 






ej2(h^) 


( 6 . 12 ) 


which vanishes as cx) due to Eq. iHD and j 2 > J 4 > 0. The same conclusion can be 

drawn for all j = n — A: with fc > 4, so that we can state that for ^ oo 




2\-l 


dvr/ 7o(n)-b (as/47r)7i(n)’ ^ ^ . (6.13) 

0, ifjVO 


Then, in accordance with Eq. (Enni), the new (non-polynomial) asymptotic distribution am¬ 
plitude is given by 






^n,0 


n>2 


7o(n) -b (Q;g47r)7i(n) 


Cy(2x - 1 ) 


. (6.14) 
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Expanding this expression in a* and retaining only the 0(Q!s)-terms, we obtain 




' 01 . 






n>2 


7o(^ 


+ 0(a^) 


(6.15) 




FIG. 9: Evolution of the pion DA for fixed as = 0.5 and three active flavors. As a nonperturbative input at 
= 0.25 GeV^ we take the asymptotic DA tpasix) in the leading order, (a) The LO asymptotic DA, (pa,s{x), 
(solid line) in comparison with the exact NLO asymptotic DA Os) (Eq. (10.1411 . dashed line), (b) The 

relative NLO corrections to yas(a^, Q^) at = 4 GeV^ are shown in units of as/(27r) following for an 
easier comparison: the exact RG prediction, Eqs. (I0.9all -- ll0.9cll . is denoted by a solid line and the result of 
0, Eq. dO.lOall . is represented by a dashed line. 


Using the explicit form of —see Eqs. —this expression can be recast in 

the form 


CXa 


= (Pas{x) 1 + ^ {CF(j)F{x)+ bo(j)bo{x)) +0(tt^) (6.16a) 

L 47r J 

which coincides with Eq. (31) derived in Ref. [^. Here, the following abbreviations have been 
introdnced 


(j)F{x) = 


n>2 
2 


4A 


' 4 (2n + 3) 

(n + l)(n + 2) [?7,(?7, + 3) 


+ Ano - ii{n + 2) + i>{l) 


Cl/\2x-l) 


= In 

Ki.^) = Y 


- - 1 

X 




n>2 


- -- -Cl^^{2x-l) = In [x(l - x)] + - . 

n(n +l)(n + 2)(n + 3) ^ ^ ^ ^ ^^3 


(6.16b) 

(6.16c) 


Hence, we conclnde that in the fixed as case onr approach fnlly reprodnces the two-loop 
resnlt of Ref. and the part (j6.16cjl coincides with Mikhailov’s finding in [i^. But our 
approach contains more information. Because of the RG improvement, Eq. (j6.14fl contains 
the coefficients Ln,o instead of Mo,n entering this equation in the approach of [2|— cf. Eq. 
fj6.15|) . The advantage of the coefficients Ln,o is that they are much smaller, as we clearly 
see from Fig. IHlby comparing the two approaches. In Fig. El we show the comparison of the 
two approaches for the pion DA itself. Note that in our analysis we use the same numerical 
values as = 0.5, Nf = ?>, jj? = 0.25 GeV^, and = 4 GeV^ as in and take into 
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account in Eq. the first 100 nontrivial terms (n = 2,4,..., 200) to approximate the 

evolved distribution amplitude. It is important to stress that our results, obtained in the MS 
scheme of perturbative QCD for a hxed coupling, can be related to those derived before by 
Muller [ 2 ^ in the conformally covariant subtraction (CS) scheme. More specifically, the 
solutions of the ERBL evolution equation in each scheme can be interlinked in the conformal 
limit (/9 = 0) on account of a finite refactorization 21| . This means that our solution (lfi.lba|l 
in the MS scheme and Muller’s corresponding results in the CS scheme 1^ are connected 
by a RG transformation matrix determined by the special conformal anomaly 0 - 


VII. CONCLUSIONS 

In this paper we analyzed the ERBL evolution equation for the meson distribution ampli¬ 
tude at the two-loop level. Our main conclusions can be summarized as follows. 

• We worked out a procedure to solve the ERBL evolution equation in the NLO approx¬ 
imation, Eqs. innnii . which ensures the RG properties. 

• Using this method, we obtained a RG-improved solution of the NLO ERBL evolution 
equation in the form of quadratures, Eqs. (HU), (I43E1) . 

• We worked out an approximate version of our procedure, given by Eqs. (ICTIll and 
(gTTH . which is completely analytical. Its accuracy is rather high with an uncertainty 
of the order of only a few %, except for the case of ^ 4 , 0 , for which the error is about 
- 10 %. 


• We analyzed the importance of the NLO evolution for the inverse moment of the meson 
distribution amplitude. We conhrmed the conclusion drawn in that it is possible 
to use the LO evolution for this quantity with the induced overall error being smaller 
than 1%. 

The most discernible theoretical result of this investigation is that the RG-improved evolution 
is slower and the higher harmonics stronger suppressed relative to the approximate MR 
scheme, enhancing the self-consistency of QCD perturbation theory. 
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APPENDIX A: BETA-FUNCTION AND ANOMALOUS DIMENSIONS 


The standard /3-function coefficients are given by 

, llN,-2Nf , 34 

bo = -^^ ; b 


10 


i = y^'c - (2UF + yiV,)iV^. 


(A.l) 


To one-loop, the anomalous dimensions read 

n+1 


7o(n) = 2Cf 


4^--3- 


2=1 


{n + l)(n + 2 ) 


(A. 2 ) 


whereas the anomalous dimensions for arbitrary n at the two-loop order may be found in ^ 

na, with 

34450 830,^, 662846 31132 ^^ ,, 

A(0) = 0, 7,(2) = 71(4) = yyy - . (A.3) 

The normalization of the Gegenbauer polynomials 'ipn{x) with respect to the weight function 
r 2 (x) = 6 X (1 — x) is 


(A.4) 


N’' = I +" 3 ) ■ 

The values of the first few matrix elements Mj^n have been calculated numerically in to 
be 


Mo 2 = -11.2 + 1.73Nf, Mo 4 = -1.41 + 0.565Aj, M 24 = -22.0 + 1.65Nf . (A.5) 


In particular, for Nf = 3, they read 


Mo 2 = -6.01, Mo 4 = 0.285, M 24 = -17.05 . (A. 6 ) 

It is worth pointing out here that Miiller in Ref. [a has obtained analytic expressions for the 
matrix elements „ for all values j = 0,2,... <n = 2,4,... ^ 


Mj,n 

^( 1 ) 


A 


n,3 



Cnj [7o(’^) - 7o(j)] , 


(2j + 3) 


- 70 (j) - 2bo + SCpAnj 2CF{Anj - tl^jn -F 2 ) - 0 ( 1 )) 

2(n - j) {n + j + 3) (j + 1) (j + 2) 


i! 


^+J+A '^ 



+ 2%lj{n 


j) -'ip{n + 2) . 


(A.7) 

(A. 8 ) 

(A.9) 


APPENDIX B: PROOF OF Eq. dOm 

Using the principle of mathematical induction we want to prove that 

^r^,r^-fc(Q^/l^) = '^n,n-k{Q^) - '^n,n-kiA^) - (Q^/l^)4^n-i,n-fc(h^) (B.l) 

0<j<k 


The function ijj{z) is defined as usual by Jp{z) = InT/z). 
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is the solution of Eq. (ITTHl) . 

[T] We proved this for k = 2. 

Let us assume that it is valid for all A; < m — 2 and then prove that it is valid for 
k = m. We have 


= —Zn^n-ra{Q‘^ 1 Zn^n-raiQ'^ iff') Zn^n-ra{(f‘1 — 


0<j<m 




0<j<m 


= ' [Zn,n-j{Q^, k^) - Zn,n-j{Q^, “ ^n,r^-i(g^ k^)] ^n-j,n-m{k^) 

0<j<m 

+ Y Zn,n-jiQ^,q‘^) [^n-j,n-miq'^) -'^n-j,n-m{k‘^)] ■ 

0<j<m 

Now we can use Eq. flTT^ to rewrite the last line in the previous equation in the following 
way (on account of j < m — 2 in all sums) 

n—j,n—m(^q ) ^n—j,n—m(^k ) ZYi—j^n—mi^q i k ') ~k ^ ^ Z-^—j ^i—ki^q ,p n—k,n—m(^k ) • 

j<k<m 

Then we hnd 

Z{m, k^) = Y [^n,n-j(Q^ k^) - Zn,n-jiQ^, 9^) “ Z^^n-M^, k^)] n-j,n-m{.k^) 

0<j<m 


+ Zn,n-j{Q^, q‘^)Zn-j,n-m{q‘^, k"^) 

0<j<m 

+ Zn^n-kiQ"^ iq^)Zn-k,n-j{q^ 1 k‘^)^n-j,n-m{k^) ■, 

0</c<m k<j<m 

where in the last line we have introduced new indices {j,k —> k,j). Changing the order of 
summations in the double sum gives 

Z{m,Q‘^,q^,k^) = Y [Zn,n-jiQ‘^,k‘^) - Zn,n-j{Q‘^,q^) - Zn,n-j{q^, k^)] n-j,n-mik‘^) 

0<j<m 

+ Zn,n-j{Q'^, q‘^)Zn-j,n-m{q'^, k^) 

0<j<m 

+ Y ^n,n-k{Q‘^,q^)Zn-k,n-j{q^,k^)'^n-j,n-m{k^) = 

0<j<m 0<k<j 

Y ' [Zn,n-j{Q‘^, k^) - Zn,n-j{Q^,q^) “ Zn,n-j{q^,k^) 


0<j<m 


+ Zn,n-k{Q‘^ ■,q^)Zn-k,n-j{q^ 1 k^) ^n-j,n-m{k^) 


0<k<j 


+ Y Zn,n-j{Q^,q^)Zn-j,n-m{q‘^,k^) ■ 

0<j<m 
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By virtue of j < m — 2 in all sums, we may use Eq. (ITTHl) for the expression inside the square 
brackets to show that it is identically equal to zero. Therefore, we find 


= Zn,n-m{Q‘^, - Zn,n-m{Q‘^,q^) “ 

= ^ Zn,n-jiQ^,q‘^)Zn-j,n-m{q‘^,l^^) ■ 

0<j<m 


This way we obtain just the desired expression, cf. Eq. (IXTHl) . 
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